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StRT II THRUST CORRECTION FOR DOUBLY CHARGED IONS 
Ralph R. Peters 


The thrust of the SERT II thruster measured in space has been shown 
to be slightly less than that calculated from electrical measurements on 

n 1 

the thruster.^ Since the electrically based thrust was corrected for 
ion beam divergence one would expect this thrust difference to be due to 
doubly charged ions. In order to verifjr this postulate the ion beam of a 
SERT II thruster was probed using a mass spectrometer suitable for double 
and single ion discrimination and the thrust error due to the double ions 
present was determined. 

Procedure 

A SERT II thruster was installed within the Colorado State University 
vacuum test facility and operated at the design 100% thrust level. The 
thruster had been modified in that it had electromagnets and independent 
main and cathode flow systems, but these were controlled in such a way 
that in-flight operating data were reproduced to the extent indicated in 
Table I. 

Table I - SERT II Operating Conditions 



In-flight 

Current Test 

Arc Current 

1 .7 A 

1.7 A 

Arc Voltage 

37 V 

37 V 

Bearn Current 

253 mA 

258 mA 

Mass Flow Rate 

313 mA 

307 mA 

Screen Voltage 

3000 V 

3000 V 

Accel Voltage 

-1500 V 

__ -1500 V 




Wlins 0|»i',itin-| at the test roniUtloiir. indicated the ion heaiii was 
prohed with an t x U momentuia analyter'-^] and the resultant data „e,-e 
analyzed in the manner described in Appendix A to obtain the double and 
single ion current contributions to the measured beam current. 


Results and Conclusion 

The double-to-single beam ion current ratio U**/l*) determined for 
the SERT II thruster at the Table I operating conditions was 7.8?„ The 
ratio of true thrust to that which would be calculated assuming only 
singly charged ions and operation at the measured beam current is given by: 


Double Ion 
Thrust Correction 
Factor 



r+ + 


1 ^/ 2 ^ 


(1 + 


( 1 ) 


Substitution of the measured double-to-single current ratio into 
this expression yields a double ion thrust correction factor of 0.98. 
The SERT II thrust level as determined from miniature electrostatic 
accelerometer measurements in space was 27.4 niN ± The thrust 

calculated from electrical measurements and corrected for ion beam di- 
vergence was 28.2 mN t 2.2S.[0 Multiplying this calculated thrust by 
the double ion correction factor yields a thrust of 27.6 mN which agrees 
with the measured thrust to well within the accuracies involved. 


DOlJlil.Y UIAKGLI) ION I’KOUUGTKiM tlM.IKON 
BOMBARUMfNT THRUSThRS 
Ralph R. I’etf'r , 

Doubly charged ions are generally considered an undesirable component 
of an electron bombardment thruster plasma because of the sputtering 
damage they cause. In order to select design and operating constraints 
which minimize doubly charged ion concentrations it is desirable to have 
an accurate model describing the mechanisms by which they are produced and 
lost. This work describes such a model and its verification through tests 
on the 15 cm diameter SERT II mercury ion thruster. Although this particu- 
lar work is applicable to mercury discharges only, the procedure should 
apply equally well to other propellants. 


Theoretical Model 

The reactions which are considered important in the production and 
loss of ions in mercury, electron bombardment ion thrusters are illustrated 
in Figure 1. The symbols used in the figure represent the following specie 

Hg - nsutral ground stato fnorcury 

Hg"^ - metastable neutral mercury (6^P and 6^P,, states) 

Hg - resonance state mercury (6^P, and 6^P, states) 

+ 

Hg - singly ionized ground state mercury 

Hg - singly ionized metastable mercury and 6^0.,., states) 

++ -5/^ b/c: 

Hg - doubly ionized ground state mercury 

The arrows in Figure 1 indicate the various interaction routes con- 
sidered in this analysis. They show for example that ground state singly 
charged ions can be produced from neutral around state, resonance state, 
and metastable state atoms and that they .an l)e lost as a result of single 
ion migrcition to the l.oundary and produ. tion of me ta stable doul)le 
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ions and doubl(> ions by olectron bombardiiienl. Fho intoniiediate rosonanre 
and nietastablo states used in the analysis were Sv^lected because they have 
significant electron impact cross sections of formation over the electron 
energy range of interest and sufficiently long effective lifetimes so that 
they can enter into ion production processes before they decay. In going 
to the plasma boundary an excited atom or an ion could be going to a dis- 
charge chamber wall where it would be neutralized and returned to the 
discharge as a neutral ground state atom or it could be going through a 
grid aperture in which case it would be extracted from the discharge region 
and replaced by a neutral atom from the propellant feed system. In either 
case this process represents a loss rate for any of the excited states. 

These losses to the boundary are indicated in Figure 1 by the dotted lines 
to the wall of the chamber and then a large arrow back to the neutral ground 
state. As suggested by Figure 1 resonance state atoms can also be lost by 
photon diffusion. These atoms have short lifetimes for spontaneous photon 
emission but the emitted photon is readily absorbed by an adjacent ground 
state atom thereby producing another resonance state atom. Since the 
transport time of the photon is i.. all compared to the excHed state life- 
time the excited state can be considered to exist continuously. The photon 
can however eventually diffuse to a boundary where it will be lost. This 
is equivalent to the loss of a resonance state atom; a loss mechanism 


represented in Figure 1 by a dotted line conveying a photon to the wall and 
a branching line going from the resonance atom to the ground state atom. 


The determination ot the various specie densities requires the 
equating of production and loss rates for eai h specie in accordance with 
tfie above model and as outlined in Reference [3]. Where possible, cross 


sections for the reactions 


implied by figure 1 W( 0 'e selei ted fi'om published 


I 


(r 

expor iiiK'iiiii I ■*’ ^ Jl I > i ir, , .!•. i iiiii'. iiot iiviril'- 

alVli', l‘.h(!Oi'(' I, i(.a 1 cioss ‘.cf.i, iuir; wiTf cilhcr olil ,i i ih’'l iiimi l.lu' I i tf>ra1;urtJ ^ ^ 

Ini 

or Cci 1 r.u I a tfid iriiii'i l.lio (iry/iir.ki iipproxinial.ioii. " All t(Mi i.ions except 
tliose occurr i ii(| at the itlasnia houii'iary wei'(> ion , KiinaKl to la' induced by 
electron boinbardinent. the l)()iiibard i rip elec Iron (iinup wa', assumed to consist 
of a Maxwellian component (defined by a density and temiierature) and a 
monoenerpetic or primary component (defined l),y a density and energy). 

The model assumes electron induced reactions occur only within the 
primary electron region of the thruster defined by the surface of revolu- 
tion of the critical field line and screen grid because of the higher 
electron energies and densities that exist in this region. Diffusion loss 
effects are on the other hand deriendenc on the conditions at the surface 
of this region. In order to obtain accurate results from this model it 
was necessary to input to the computer program aptiropriately weighted 
volume and surface area averaged eleett'on propcrrl ies . Tire equations used 
in the model are described in detail in RelorcnK.e [3J. 

Procedu res jnid_ Resuj ts 

Verification of this theoretical iiioilel of double; ion production and 
loss mechanisms re(|uires essentially simultaneous mea'.urement of the plasma 
properties within an operating thruster and the siiigli' and double ion con- 
tent of tile ion beam. These iiieasut'eiiients wen; made nti a fd.RT II thruster 
operating at three different arc volt.ige loiulitioir. .md witli two different 
sets of grids (STUT II and higli porveance dished ai'ids). Ihc plasma proper- 
ties were iiieasms'd usitni the movable I angmuir pinin' and recording system 
desirilH'd in K('t('ren<(' 1 '> | . llie a . ,oi iated 1 angmuir probe traces wore 
rnlidited V/ i th lire prohe at si’Iceii dilleriii) pninis in 1 he di'iihai'ge rhaiilbet 
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and the data were then analyzed usinq tlm 
tleference [10]. 


nuiiK'r ical method doscrihed iti 


An E X B momentum analyzer 
single ion currents in the beam 


was used to measure the total double and 
in accordaine with the procedure for data 


acquisition and analysis described in Appendix A. The integrated double- 
to-stngle ion current ratios obtained using this procedure are listed under 
the measured operating data section in Table I. 

A typical set of plasma property profiles obtained with the standard 
SERT II grids and during operation at a 37V arc voltage is presented as 
Figure 2. Also identified is the critical field line determined from iron 
‘ngs maps for this thruster. Volume and surface area averaging of these 
properties over the primary electron region defined by this field line 
yields the values designated as computer program input parameters in the 
second row of Table I. These measurement, analysis and averaging pro- 
cedures were repeated for each grid and operating condition investigated. 
The resultant averaged plasma properties, which are listed in Table I 
together with the geometrical properties of the primary electron region 
and the grids were then used as input to the computer model which has been 
described. The densities of the ionic and excited atomic species considered 
in the analysis were computed, these densities were normalized with the 
total heavy particle density in the thruster and they are listed as cal- 
culated normalized densities in Table I. The results of Table I show for 
example that the thruster operating with flat grids at 37V arc voltage 
would be predicted to have 68/. neutral ground state atoms, lit: neutral 
resonance state atoms, 6.8/ singly charged ground state ions and 0.2 
doubly charged ground state ions. 
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I iiid lly Uio C(i |(:iilrt1;ocl production rates for s iii'ily and doubly charged 
ions produced throuf)h tlie various inteniKMiiate states and tionnalized by 
the total production rate lor th(> specie are indicated in Tab'o I along 
with the fraction of the associated interactions effected by primary elec- 
tron- (in parenthesis), for example at the 37V, flat grid operating point 
60;!'. of the single ions were produced as a result of electron interaction 
with neutral ground state atoms. These neutral ground state-to-single 
ionic interactions were induced by primary electrons 22% of the time and 
by Maxwellian electrons the remainder (78%) of the time. 

The data of Table I show the resonance and metastable atom states 
are important intermediate states for the production of single ions. They 
also show that while the singly charged ionic ground state is the most 
important intermediate state in the double ion production chain, its im- 
portance drops off at higher arc voltages where electron temperatures and 
energies become sufficient to ionize ground state neutral mercury directly 
The measured double-to-single ion current ratio given in Table I 
divided by the guantity (2 /2) is equal to the double-to-single ion density 
ratio prevailing in the discharge region. The value of this density ratio 
is compared to that calculated using plasma properties in Figures 3 and 4 
for the operating configurations identified in Table I. The double-to- 
single ion density has been plotted against both propellant utilization 
and arc voltage in these figures and the figures sliow the most consistent 
correlation occurs with propellant utilization. The agreement between 
tlieory and experiiiumt is observed to be within 

because th«' dislu'd grid thruster operates at high(>r flow rates and 
beam i.urrenls thi> al)solute iiia.in 1 *ud(’ ol the doul>1e ion density is also 
. i()ti i r i c<int I y hii|li('r in tliis tlnuslet. 1 lie extent ol th(! difference in 
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double ion density predicted by the computer model for the two thruster 
grid systems is shown in Figure 5 and is seen to be 2 to 5 times higher 
in the dished grid thruster. 

conclusion 

The mercury discharge model and analysis presented yield double-to- 
single ion density ratios which agree with experimental measurements to 
within 25%. The analysis shows that doubly charged mercury ions are 
produced predominantly via the singly ionized ground state with direct 
ground state neutral -to-double ion production becoming more significant 
at high propellant utilizations. Singly charged ions are produced in 
significant numbers from intermediate metastable and resonance states in 
addition to direct ionization from ground state neutrals. 




SINGLE CUSP MAGNETIC FIELD THRUSTER 
John R. Beattie 

The single cusp magnetic field (CMF) discharge chamber geometry has 
been investigated as a viable approach to achieving a uniform current 
density in the exhaust beam of an electron-bombardment ion thruster. The 
importance of achieving a uniform profile and the single cusp magnetic 
field concept are discussed in detail in References [11 and 12]. 

The cusped magnetic field geometry was chosen on the basis of applying 
a neutral residence time criterion to a divergent magnetic field geometry. 
This approach suggested the beam profile of a divergent field chamber could 
be made more uniform by increasing the ionization probability at the thruster 
periphery and by confining the axial magnetic field to the chamber wall. 

This would result in a uniform ion density in the vicinity of the screen 
grid and, as a result, a uniform ion beam profile. 

The effectiveness of the uniform ion density approach to achieving a 
uniform beam profile has been demonstrated in the past by the radial field 
thruster*^^^^ and more recently by the multipole thruster. The beam 
profile flatness parameter F (defined as the ratio of average to maximum 
beam current density) of the radial field design was 0.67 using mercury as 

the propellant while the multipole design was in the 0.70 to 0.75 range 

using argon and xenon as the propellants. Typical values for the SERT II 

divergent field thruster are 0.40 to 0.49. Faraday probe data presented 

in this report will show the cusped field discharge chamber has a beam 
flatness parameter in the 0.70 to 0.82 range. 

Ap paratus 

A SERT II discharge chamber was modified to produce the cusped 
magnetic field geometry illustrated in Figure 6* The (hanges consisted 







ciisPKu MAc.NrTic ririj) thruster 
SCHEMATIC 
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of the? addition of a centor iiia()m>t pole? piece? and thf» installation of a 
movable rear anode in the upstream end of the discharge chamber. The 
thruster length-to-diameter ratio (L/D) could be varied by replacing the 
downstream section of the chamber. Additional modifications included the 
replacement of the SERT II permanent magnets with independent front and rear 
electromagnets and the installation of a variable magnetic baffle. By re- 
versing the polarity of the magnet power supplies, the cusped magnetic 
field geometry could be converted to a divergent field geometry. Iron 
filings maps of the cusped and divergent fields are presented in Figure 7 
for a chamber length-to-diameter ratio of 0.53.* The SERT II accelerator 
system was replaced with a compensated dished grid optical system and the 
propellant feed system was modified to permit independent control of the 
main and cathode flow rates. The main and neutralizer cathodes were re- 
placed with 0.5 mm diameter orifice cathodes to facilitate sustained 
thruster operation at high beam currents. 

A movable Langmuir probe was installed inside the discharge chamber 

to obtain plasma diagnostic infonnation. In addition, a Faraday cup probe 

which could be swept through the ion beam at any desired axial location 

was installed to obtain beam current density profiles. An ^ x momentum 

[ 2 ] 

analyzer was installed in the vacuum facility to measure the charge-to- 
inass composition of the ion beam. The analyzer installation was designed 
to permit the instrument to be swept throuqli the ion beam with the probe 


Tile length-to-diai'ieter ratios relerred to in ttiis *'('port arc O.iM, 0.30, 
and 0.53. These differ sli(|titly from th(» previously reporftnl valiu's of 
0.3, 0.4, and 0.6. The differetui' is due to a slightly difterent di'fini- 
tion of tlie Icngtii L and does not imply a gi'oiiietrical d i f I (M'em (>. 
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pitclu?(l tit: various atuili's, llris (icniii i.l.oil I, tin siii(| |o and doiihio ion dis- 
porsion protilns to hi' dotoniiinod as a luiicl.iori of ttirustor radius. 

All data iirnsonlu’d vans' olitiiiiind at i |kV and “O.AkV soroi'ii and 
accolerator tjrid |iotont i.i Is and O.tA toopor currc’nt. Solid symbols tiavo 
been used on per foriiianco curves to indicate the' 37V arc voltape condition. 
Discharge power cal cul atioiv> ir ludi' the keeper losses and the total propel- 
lant flow rate was 730 niA unless indicated otherwise. 

Langmuir probe data presented herein were analyzed by the numerical 
procedure of Reference [10]. ^ x ll momentum analyzer data were recorded 

and analyzed using the procedure outlined in Appendix A. Details of the 
Langmuir and Faraday probe installations are given in References [9 and 15]. 
A schematic diagram of the ^ x 5 momentum analyzer circuit is presented in 
Reference [14]. 

Beam Profile 

The beam profile flatness parameter measured at a distance of 6 mm 
downstream of the accelerator system is presented as a function of chamber 
length-to-di aineter ratio in figure fl. Beam profile data reported in 
Reference [11] which were obtained at a distance of 50 mm downstream of the 
grids are also presented. The trends are seen to be the same for both 
probe locations and the two curves are shifted Ijy almost a constant amount. 
The data recorded close to tlie aaelorator grici indiiate a 10 increase in 
the beam flatne'c, paramo li'r was coali.s'd with Ifu' '.liorfi'r di'. charge chambr-r. 

Ion lieam current d(.'tisity protih's moa'anr.'d i< mm downstream of tlu' 
acceleratoi griil ari' preconted in I itiure u lor llie 'd I^T 11 and i iisped 
field tiiru'.teo'; . Ihi'M' data wore uonna I i /i-d ■ in li that tho intogratod 
Ih.'iIiii current is tlie '.amo lo> o.k h ptoliP'. Wla'M nm niii i i ,'oil in this iii.iiiuot , 
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the reciprocil of the peak normalized density Is numerically equal to the 
flatness parameter. The data of Figure 9 indicate a significant improve- 
ment in beam profile flatness parameter has been realized with the cusped 
field discharge chamber design. The beam flatness parameter for the cusped 
field thruster having a length-to-diameter ratio of 0.23 is seen to be 40% 
greater than the SERT II value at the same operating conditions. 

The effect of magnetic field geometry on beam profile is illustrated 
in Figure 10. These data are presented for the cusped and divergent magnetic 
field geometries with a discharge chamber length-to-diameter ratio of 0.23. 
The effectiveness of the cusped field design in achieving a more uniform 
beam current density is evident when comparing these profiles. The diver- 
gent field data of Figure 10 and the SERT II data of Figure 9 indicate the 
short chamber divergent field geometry has a flatness parameter which is 
27% greater than the SERT II value. Beam profile comparisons for the cusped 
and divergent field geometries were found to have the same trend as shown 
in Figure lOfor the intermediate and long chamber lengths. The effect was 
less pronounced for the long chamber (100% of SERT II length) as one might 
expect. This occurs because the magnetic field is largely axial over the 
entire chamber cross section for both configurations as illustrated in 
Figure 7. As the cusped field chamber length is reduced however the axial 
component of the field is significant only near the chamber wall. This 
results in a uniform plasma density over much of the downstream region of 
the discharge chamber. 

The axial location of the movable rear anode was found to have little 
effect on the beam profiles of either the cusped or divergent field 
geometries for all iength-to-diameter ratios investigated. This is to be 
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expected r>ince the anode position primarily determines the chamber volume- 
to-area ratio and had little effect on the plasma density profile at the 
screen grid. 

The effect of throttling on the beam current density profile is seen 
by comparing Figures 10 and 11 for the cusped field geometry having a length- 
to-diameter ratio of 0.23. The improvement in the beam flatness at the 
throttled flow condition is thought to be due to the presence o^" fewer 
double ions at the lower ion beam current. Since double ions are produced 
mostly from single ions, the double ions are preferentially produced near 
the center of the discharge chamber where the plasma density is highest. 

This condition combined with the larger Bohm velocity near the chamber wall 
suggests the double ion density profile would be more peaked than the single 
ion profile. Experimental measurements of the double and single ion current 
density profiles indicate that this is in fact the case. The argument pre- 
sented above would also explain the dependence of the cusped field thruster 
beam flatness on propellant utilization at a fixed propellant flow rate. 

The effects of throttling and propellant utilization on double ion content 
have not been studied experimentally to verify this theory however. 

Performance 

The cusped magnetic field thruster was operated initially using an 
uncompensated dished grid accelerator system having the dimensions listed 
in Table II. Performance data for this thruster configuration were reported 
in Reference [H] are reproduced in Figure 12 for comparison with data 
obtained during tlie current reportin<| period. Figure 13 presents pertormance 
data obtained with a set of compensated dished (|rids havin() the dimetv ions 
listed in Table 11. Botli Fi<|ures 12 and 13 indicate a substantial improvement 
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TABLE II. ACCELERATOR SYSTEM DIMENSIONS 



Uncompensated 

Compensated 

Screen 

Accelerator 

Screen 

Accelerator 

Thickness 

0.46 mm 

0.46 mm 

0.38 mm 

0.51 mm 

Hole Spacing 

2.5 mm 

2.5 mm 

2.08 mm 

2.08 mm 

Hole Diameter 

2.0 mm 

2.0 mm 

1.91 mm 

1 .6 mm 

Open Area 

1 

58% 

58% 

76% 

54% 


in thruster performance was realized as a result of reductions in the dis- 
charge chamber length. A comparison of Figures 12 and 13 reveals signifi- 
cant reductions in baseline discharge losses were realized with the compen- 
sated grid set which has a thicker accelerator grid with smaller hole 
diameters and the larger open area screen grid. Also, the performance data 
for the chamber having a length-to-diametei' ratio of 0.30 indicates a marked 
difference in the shape of the curves at high propellant utilization. The 
reason for this apparent discrepancy is not known. Since the data of Figure 13 
exhibit a more consistent trend with discharge chamber length reductions, the 
data of Figure 12 are felt to be the source of this discrepancy. The apparent 
error in Figure 12 could be due to an error in the flow rate data for the 
chamber length-to-diameter ratio of 0.30. If the measured flow rate was high 
the calculated propellant utilization would be less than the true value and 
the performance curve would be shifted to lower utilizations. This would 
explain both the presence of the "knee" in the data of Figure 12 and the dif- 
ference in baseline losses for the short and intermediate chamber lengths. 

The effect of discharge chamber length on performance for the diver- 
gent magnetic field geometry is illustrated in Figure 14. The baseline 
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discharge losses are seen to decrease with a reduction in chamber length 
as was observed with the cusped magnetic field geometry. The increase in 
propellant utilization observed with the shorter discharge chamber appears 
anomalous at first glance, bi*- is in agreement with both neutral loss rate^^^^ 
and double ion production^^^ theories when the volume-to-surface area ratios 
of the three configurations are taken into account. Iron filings maps of 
the magnetic field existing in the three configurations indicate the primary 
electron region volume-to-surface area ratio increases as the chamber length 
is reduced. Since theory predicts an increase in double-to-single ion density 
ratio and a decrease in neutral loss rate with increasing volume- to-area 
ratio, the shifting of the performance curve to higher propellant utiliza- 
tions is expected. The volume-to-area ratio of the short and intermediate 
chambers was found to be about the same and the crossover of the two curves 
is probably due to differences in the plasma properties existing in the two 
chambers. 

Perhaps the most notable difference between the cusped field and 
divergent field performance data is the presence of a "knee" in the diver- 
gent field data and the lack of one in the cusped field data. This is 
illustrated in Figure 15 which is a comparison of the cusped and divergent 
field performance for a chamber length-to-diameter ratio of 0.53. The 
cusped field data suggest a significant double ion population exists in 
this chamber at the high utilization conditions. The flattening of the 
performance curve illustrated in Figure 15 for the long discharge chamber 
is also apparent for the shorter chambers as can be seen by comparison of 
Figures 13 and 14. 

The effect of rear anode position on thruster perfoniian(e was found 
to be the same for all three discharge i liaiiiber lengtiis and both iiiagnetii. 
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field geometries. In each case the downstream movement of the rear anode 
resulted in a decrease in beam current for the same power input. This 
effect has the tendency to shift performance data to lower propellant 
utilization and higher discharge losses. The reason for the performance 
degradation is due to a reduction in both the primary electron energy and 
Maxwellian electron temperature as will be shown in the section on plasma 
properties. 

Double Ionization 

Double and single ion beam currents were obtained by the use of an 
articulating t x t momentum analyzer. Details of the data acquisition 
and reduction procedures are presented in Appendix A. The double-to- 
single ion current ratios presented in this section were measured at the 
thruster conditions indicated by solid symbols on the corresponding per- 
formance curves. Figure 16 presents the double-to-s ingle ion beam current 
ratio as a function of chamber length-to-diameter ratio for the cusped 
magnetic field geometry. Since the double- to-single ion density ratio is 
equal to 1/(2 /?) times this current ratio these data suggest a chamber length 
exists which results in a minimum double-to-single ion density ratio. Double 
ion production theory^^^ predicts a linear dependence of the double-to-single 
ion density ratio on the chamber volume- to-surface area ratio. The non- 
linear nature of Figure 16 suggests the plasma properties also vary with 
chamber length. The variation of the plasma properties with discharge 
chamber geometry and magnetic field configuration will be discussed in the 
next secL ion. 

Figure 17 presents the double-to-single ion beam current ratio as a 
function of chamber length-to-diameter ratio for the divergent magnetic 
field geometry. Here the effect of ctiamber length on the double ion 
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population agrees fairly well with the theoretical dependence on volume- 
to-area ratio which generally decreases with chamber length. The nonlinear 
nature of the data of Figure 17 again suggests however a slight variation 
in the plasma properties with discharge chamber length. 

The effect of the rear anode location on the double-to-single ion 
density ratio is seen from Figures 16 and 17 to be essentially the same for 
each chamber length and magnetic field geometry. The reduction in the 
double-to-single ion density ratio associated with the downstream anode 
position is consistent with the reduction in the volume-to-surface area 
ratio caused by the anode movement. However, the effect illustrated in 
Figures 16 and 17 is not due entirely to volume-to-surface area ratio 
changes. Langmuir probe data presented in the next section indicate the 
downstream movement of the rear anode had a large effect on the Maxwellian 
electron temperature and primary electron energy. 


Plasma Properties 

The cusped field thruster discharge chamber theory^^^’ has 
suggested electrons have the tendency to be trapped in a magnetic bottle 
region which exists in the upstream end of the chamber. This theory was 
verified for the three different discharge chamber lengths by Langmuir 
probe measurements. The variation of the plasma electron density with 
axial distance is presented in Figure 18 for each discharge chamber length 
and both the cusped and divergent magnetic field geometries. The parameter 
r/r^node these curves is the nondiinensional radius at which the 

probe sweep was made. Comparison of the cusped and divergent field data 
clearly indicates the axial density profile of the divergent field geometry 
is quite uniform while the cusped field profile is highly nonuniforni. The 


ELECTRON DENSITY X 10"'® (cm"®) 


- 36 - 


50, 


30 


OL 


iCENTER POLE PIECE 

k «CUSPED MAGNETIC FIELD 

a DIVERGENT MAGNETIC FIELD 

•'''’anode ® 

'■/^anode = 0*88 

A. L/D =0.23 



50r 


30 



*^^*"anode * 0*59 
'^''•^anode * ^*88 



C. L/D = 0.53 


PLASMA PROPERTY VARIATIONS 


FIGURE 18 


-37- 


cusped field profiles indicate the presence of a localized region where the 
electron density is as much as ten times the value near the screen grid. 

The axial location of this localized region of high electron density cor- 
responds to the position of the center magnet pole piece. Figure 18C indi- 
cates an increase in the plasma electron density near the screen grid for 
the divergent magnetic field geometry. This increase occurs because the 
probe crosses the critical field line and enters into the primary electron 
region as it moves downstream. This can be visualized by considering the 
iron filings map of Figure 7B. The Maxwellian electron temperature and 
primary electron energy were found to exhibit trends similar to those 
illustrated in Figure 18. 

Plasma property surfaces are presented in Figure 19 for the cusped 
and divergent magnetic field geometries with a discharge chamber length- 
to-diameter ratio of 0.23. Comparison of the cusped field and divergent 
field surfaces indicates the cusped field geometry generally results in 
much more uniform plasma properties than the divergent field geometry. 

The reduction in the radial plasma density gradient is clearly evident in 
the cusped field design and this is reflected in the more uniform beam 
profile of this thruster. The magnetic bottle effect is apparent when con- 
sidering the plasma property variations near the center magnet pole piece 
of the cusped field geometry. The general trends evident in Figure 19 
were also observed when comparing the surfaces for the intermediate and 
long discharge chambers. 

The effect of the rear anode position on Maxwellian electron tempera- 
ture and primary electron energy is presented in Figure 20 for the cusped 
field geometry with a length-to-diameter ratio of 0.23, These data indicate 
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a large reduction in the electron temperature occurred throughout the dis- 
charge chamber as a result of moving the rear anode to the downstream 
location. The reduction in electron temperature is thought to be due to 
a depletion of the high energy (large cyclotron radius) electrons in the 
Maxwellian tail as the anode is moved close to the magnetic field lines in 
the upstream end of the discharge chamber. The downstream movement of the 
anode also caused a reduction in the primary electron energy level. The 
primary energy effect appears to be more localized than the temperature 
effect probably due to the longer path length of the primary electrons for 
Coulomb collisions. The effect of anode position on the electron temperature 
and energy as described above was found to be the same for all chamber 
length-to-diameter ratios and both magnetic field geometries. 

Langmuir probe data were used to calculate the volume averaged plasma 
properties in the manner suggested in the Doubly Charged Ion Production 
section of this report The results are presented in Table III for all 
three chamber lengths, both magnetic field geometries, and both rear anode 
positions. The average plasma properties can be used to explain most of 
the double-to-single ion density ratio trends of Figures 16 and 17. For 
instance, the Maxwellian electron temperature, primary electron energy, and 
volume-to-surface area ratio are almost invariably decreased as the anode 
is moved to the downstream position. This should result in a reduction in 
the double-to-single ion density ratio and the experimental results of 
Figures 16and 17 confirm this for all cases except the short chamber cusped 
field configuration. The variation of the experimentally determined double- 
to-single ion density ratio with discharge chamber length and magnetic field 
geometry is not as easily correlated with the average quantities presented in 
Table III. Those cases which seem inconsistent are probably due to differences 
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TABLE III. VOLUME AVERAGED PLASMA PROPERTIES 




Cusped Magnetic Field 

Divergent Magnetic Field 


Anode Position 

Anode Position 

L/D ; 

* 

Property 

Upstream 

Downstream 

Upstream 

Downstream 

0.23 i 

(eV) 

4.7 

3,05 

4.9 

3.8 

I" ! 

5p (eV) 

29.4 

28.0 

29.7 

27.6 

1 

npXl0‘^°(cm“^) 

0.94 

0.76 

1.03 

0.92 

(i 

n^xl0“''°(cm“^) 

26.8 

28.1 

26.8 

31.7 

II 

V/A (cm) 

1.65 

1.35 

1.51 

1.0 

0.30 

(eV) 

4.8 

3.3 

5.6 

3.9 

i 

1 

! " 


28.7 

27.1 

28.9 

26.7 

! II 

1 

npXl0‘^°(cm'^) 

0.62 

0.77 

0.83 

0.99 

1 

1 “ 

n^xl0"''°(cm‘^) 

26.2 

32.7 

24.3 

29.4 

II 

V/A (cm) 

1 1.82 

1 

1 .63 

1.48 

1.36 

0.53 

T„ (eV) 

3.1 

3.1 

3.5 

3.4 

II 

f.p (eV) 

' 22.5 

24.1 

27.0 

19.8 

II 

npx10''°(cm'^) 

: 0.56 

I 

0.47 

0.57 

0.58 

! i( 

n^xlO*’“(cm‘^) 

1 28.6 

27.2 

27.7 

33.5 

It 

V/A (cm) 

1 1 .93 

J 

1.93 

1.14 

0.81 


* T - Maxwellian Electron Temperature 
tn 

■, - Primary Electron Energy 

P 

Mp - Primary Electron Density 
- Maxwellian Electron Density 

V/A - Volume-to-Surtace Area Ratio of Primary Electron Rc'()ion 
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in the plasma density existing near the screen grid. The data of Table III 
also show a general tendency for primary electron energies and densities 
to increase as the length-to-diameter ratio is decreased. 

Double Ion Correlation 

Double-to-single ion current ratios have not been measured at high 
propellant utilization due to thermal expansion and subsequent shorting 
of the grids which occurs when the thruster is operated for extended periods 
at high power conditions. The double-to-single ion current ratios which 
have been measured were obtained at a propellant utilization of about 85% 
and were found to be less than 8%. Performance measurements at higher 
power conditions suggest the ratio may go as high as 30% at a relatively 
modest discharge power loss of 600 eV/ion. This large apparent double ion 
ratio observed at the higher power conditions is observed with the cusped 
magnetic field geometry but not the divergent field geometry (see Figure 15). 

Some justification of the assumption that the cusped field thruster 
produces a large fraction of double ions at the higher power and beam current 
conditions can be made on the following basis. The shifts in the performance 
curve observed with the typical cusped field thruster as propellant flow 
rate is increased is illustrated in Figure 21. The performance curve at 
low flow rate resembles that of a properly designed discharge chamber. That 
is, the discharge losses are fairly constant until the "knee" of the curve 
is reached whereupon the losses rise suddenly and a maximum propellant 
utilization is reached at a high power condition. Increasing the propellant 
flow rate causes the performance curve to shift toward a higher maximum 
propellant utilization. As the flow rate is imreased even further the 
steep portion of the curve l)eqins to flatten out and pt iiellant utilization'. 
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from a [)erforiiiance rurvo olit.n’iKHl at a low |.iropol 1 ant flow rate where the 
double ion density is assumed to lu> rioql inihle. Then assuming a constant 
neutral loss rate at this discharge power in accordance with Kaufman's 
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model , the inaximum propellant utilization can be calculated for each 
flow rate by the expression 

"max ' ’ - "'losT'’' (2) 

where n is the propellant utilization, iIi the total propellant flow rate, 
"^loss neutral loss rate. From this maximum propellant utilization 
and the measured utilization, one can calculate a doub le-to-single ion 
density ratio assuming the difference in utilizations is due to the presence 
of double ions. The double-to-single ion density ratio can be expressed in 
terms of the measured and maximum utilizations by the expi'ession 
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The beam current can be expressed in terms of the ion densities by the 
expression 


an +2/2 n^^ 
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For double-to-singte ion density ratios less than about 12:;; the quantities 

n* t 2n*^ and n" t 2/1 n*" differ by less than about 6;;: and to a pood 
approximation one can write 
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That IS, the double-to-sing1e ion density ratio should vary almost linearly 

with beam current at a given discharge power (eV/ion) as one varies mass 
flow rate. 


Double-to-single ion density ratios were calculated from performance 
data obtained at various flow rates by tne use of L()uations 2 and 3. The 
calculations were made for discharge be.- length-to-diamete.- ,-atios ot 
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at these results were valid and tfiat the excessive propellant utilizations 
measured at the hiqh beam currents ar'e diie to tin? existence of significant 
double ion populations in the cusped fie'd discliarqe chamber. 

Conclusions 

Both the beam profile flatness and thruster |)erformance are improved 
as the discharge chamber of the single cusp magnetic field thruster is re- 
duced to 63% of the SERT II thruster length. The beam profile flatness 
parameter of this short discharge chamber is 0.70 when operating at 85% 
propellant utilization (730 mA total flow) and a discharge loss of 250 eV/ion. 
This represents a 40% increase in beam flatness over the SERT II thruster 
at approximately the same propellant utilization and discharge power level 
(eV/ion). Throttling the thruster to a 460 niA flow rate results in an 
increase in the beam profile flatness parameter to a value of 0.82. 

The measured double-to-single ion density ratio for the short chamber 
cusped magnetic field geometry is 2.4' which is slightly less than the 
value obtained for the SERT II thruster at comparable operating conditions. 
The double-to-single ion density ratio is thought to increase to much higher 
values at higher power conditions, but measur'ements to support this theory 
have not been made at the present time. 

The divergent magnetic field geometry w-r, lound to have a less ui'form 
beam profile than the cuspt'd tiffld i|eoiH('i,ry lor all tlirr’e > hamh(?r lengths. 
However, the lieaiii (trnf ile ol llu’ divergi'ol I ii'M (leometry v/as ( onsidei'ahly 
more uniform than the SI.Rf II UiriistiT piolili'. 

Tfie nos i l ion of the rear anndi' w.r. Ininul In liavi' a igniliiant cMe<t 
on the Maxwellian (.'I'm trot. I einpina t m I' . Movriniuif (O llu' risir aiindc in the 
downstream dirction n'ailird in i igniii' ml i I'llij! l inn in the dnub I ' I n- 




ION Ol’TJCS STUDY 
Orvioiiit; Aston 

An accurate knowledge of the ion beam half angle ('‘^ 1 ^^) is of con- 
siderable importance when applying ion propulsion devices to actual 
mission situations. This angle, which defines the cone confining the 
expelled ions, should be small to maxiiiii/e the thrust at a given beam 
current. It would be advantageous to know, not only what accelerator 
geometries produce the lowest half angles, but also how the half angle 
could be expected to vary as the perveance of a grid system varied due to 
power variations during a mission. 

Computer solutions are available, which predict ion optical trends 

as a function of grid and accelerating voltage parameters, for a single 

ri71 

hole geometry.'- This experiinental ion optics study was instigated to 
verify the accuracy of these solutions by comparing them with results 
obtained from actual single hole grid sets. Performance expected with 
typical multiaperture grids was then investigated by using a nineteen 
hole hexagonal array of approximately one cm^ area. This was assumed 
adequate to model the adjacent hole interactions found in full size grid 
systems. 

Parameter Definitions 

Throughout this report the geometrical grid parameters and symbols 
used are consistent with those employed liy Kaufman^^^^ and illustrated in 
Figure 23. 

These are listed below: 

-■ separation between m hmu) and a^a (’lerator grids 
d^ diameter ol aMebo alor api'i t.in c (■. ) 
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SCREEN ACC 

GRID 



t - thickiu's'. of <i(:(.(>|(.>r.ifor (|rirl 
t, ‘ tlvickiK.'ss t)f scroon <)i'id 

Kriufnian offV^ctivo atu.cloration lonqth 
defined by: 



“max ~ '’'^ximuiii ion beam half angle. 

Throughout this section perveance calculations were normalized by 

V. .. 

multiplying with the ratio (-j--)’ and grid parameters were nond imensional ized 
with the screen aperture diameter (d^). 

Apparatus 

A simple, mildly divergent field, 8 cm electron bombardment ion source 
was constructed and operated on argon propellant for these tests. Tungsten 
wire filaments were used as both the main and neutralizer cathode emitters. 
The magnetic field was derived from a long solenoidal winding extending the 
length of the discharge chamber, with a single Helmholtz coil winding 
positioned at the chamber's rear. The geometry was such that the field at 
the front of the discharge chamber was 60':. that of the rear. A cylindrical 
anode was employed and stainless steel construction was used throughout. 

The screen and accelerator’ (p’ids used were labricated from thin slieet 
graphite. This material wa'. easy to iiiadiiiie and Wti'. available as a stock 
item in the dt!sired thicknesses, live centimeter s<iuari' platr's ( on1.<ri n i m| 
the grid palter’ir. wi're positioned on a masking plate whirh i overed most 
f)t ttie down'' Lreaiii ('tul ot the discharge chamiu'r. Ilie c a |i u l.ifc'cl Viirialion 
in ion density acro'.'. the arisi ol tlio'.e grid patti'rir, w<r. |e'.' Ili.oi ■< . 
Vari(d)le '.epaisif ion o! the 'iricf, was a(|iii'Ved I'V using Inin mi' a w.i .Ikm s 
wti i ( h were rep I a i ed a I t er c sn li i la I a inn * i ■ . i V" i ' I i ' 1 1 '> I i i ' a I I a i 'a t dc iwn and 






lar<|e leakaqo ciin'eiits. 

All tests were conducted in a 30 cm diameter pyrex bell jar pumped 
by a meclianical pump in series with a 10 cm diameter oil diffusion pump. 
The power supplies and associated volt meters and ammeters were connected 
to the source in the manner shown on Figure 24A. The argon flow rate into 
the source was sensed by a Hastings flow meter and displayed on a digital 
readout. 

Ion current density profiles were measured in the ion beam using a 
movable probe rake composed of twenty 0.70 cm x 1.43 cm stainless steel 
Faraday probe sensors whose output was sensed by the system shown in 
Figure 24B. The probe rake position relative to the thruster and the de- 
tails of its construction are illustrated in Figure 25. Once the source 
had reached stable operating conditions, the rake was moved through the 
ion beam and brought to rest at that location where the current sensed 
by one of the central probes peaked out. Current to each of the twenty 
probes was then measured to determine the ion beam profile. 

Operating Conditions 

In order that as few variables as possible entered into the final 
results, the ion source was operated at the following conditions for all 
tests: 


Arc Voltage 

40 V 


Magnet Current 
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was assumed to be 5 volts. Beam current was controlled for the tests by 
adjusting the cathode current and hence its emission. From this beam current 
and the acceleration voltage and grid geometrical factors the normalized 
perveance was calculated. Through the course of the data collection the 
average bell jar pressure was in tir, low 10“' ton' range for the single 
aperture grid tests and the middle to high lO"'’ torr range with the multiple 
aperture ones . 

I on Beam M e a s u rem en t s 

Ion beam profiles measured with the Faraday iiroho rake show two ion 
populations are present in the beam--high velocily beam ions and charge 
exchange or background ions. These two populations overlap |iroduf ing an 
ion beam profile like’ ttu’ one sliown in Fi(|ure ;'P/\ and it is ik'i. essary to 
separate tlieiii so the' lieam ion lisiji'i toeie'. ran he del ined. Oin e llie true 
beam ion prol ill' Inr. Iieen idenlli ied, I in’ hall angle id the lie. mi (.an lie 
de ti.'l’liri nod . In (.lie pri".enl i .1 .r iln' iih|i( I vei'. In M’|ni t iliis hall aiigli’ 
in sui.ii a wav ttia t I 1 1 li ■ I i iii ■( I a > "n 1 . 1 1 ir . i l um, ui 1 ' j 1 na I 1 no a I I In ■ gi i '.i' 
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Separation or l.l.o ,,,,1 , Impm.' Im,,,,. io„ . Mrr,.„|., was arxoi.ip- 
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ion current density was (,,i |,.u lairsi a-, i,wi<x Ihe prodin L of the ion charge, 

Bohm velocity and (Juircie exc.iiaiKje lou (|e||^ j (,y determined from Kaufman's 
C*l 9] 

model. The factor oi tv/o was mM.ipoialed into this calculation be- 

cause it yielded a current density equal to the measured charge exchange 
current density at locations far from the ion beam centerline. The beam 
ion envelope line was then extended in the manner suggested by Figure 26B 
to intersect t^e ^harge exchange baseline and define points at the probe 
rake location which were used to calculate the desired half angle. As a 
check on this procedure several current density profiles were integrated 
along the charge exchange base! ..c out to these intersection points and 
this yielded a current that was 9(T. ot the measured beam current. This 
technique for half angle determi nation was applied to all of the multi- 
aperture data. Consistent re mj 1 1 s wen^ rT.tained over the middle to high 
perveance t ange of operation, liowevcr a i. i hi’ lower perveiuu e conditions 
some scatter was evident. In thi on,.ra, ,,a, reuime a linear extrapolation 
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poor focusing of the accelerator svM..«, „i (he low pcrvoai.ces . As a result 
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conic frustrum was shown to enclose 90% of the beam current consistently. 

The diameter of the ion beamlets emanating from a grid pair was also 
investigated. This was accomplished with a given single aperture grid pair 
by increasing the beam current until a sudden increase in accelerator impinge- 
ment current was observed. The onset of this increase corresponded to the 
beam current and hence perveance at which the beamlet diameter was equal to 
the accelerator aperture diameter. The test was repeated with different 
accelerator aperture diameters to obtain the variation in beamlet diameter 
with perveance. 

The increase in impingement current with beam current at the onset of 
direct impingement was less pronounced with the multiple aperture grid sets 
than it was with the single hole geometry; making a beam current determina- 
tion at the transition point more arbitrary. To circumvent this difficulty 
and produce consistent results, the impingement current was recorded over 
the range of beam current (perveance) variation. Plotting impingement 
current against perveance and extrapolating the straight line portion of 
the curve, a baseline impingement current was obtained. The start of direct 
impingement was said to occur at the beam current where the actual curve 
had departed from the baseline by an amount equal to fifty percent of the 
baseline current. The beamlet diameters were therefore equated to the 
accelerator hole diameter at this perveance condition. 

Resul ts 

The variation in ion beam half angle as a function of normalized 
perveance for the single aperture grids is presented in Figure 27 with non- 
dimensional ized grid separation distance as a parameter. Symbols used nn 
this figure are defined in Figure 23. Normalized perveance defined as 
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determined by Chiid's law has a maxiamni value of 
3.03 X 10-' amp/volt V. ^ 

the ratio of net-to-total accelerating voltage. Open symbols shown on 
the figure correspond to experiu.nta1 results and they are bracketed by 
etror bars to indicate the accuracy of the results at various perveance 
levels. The numerical results obtained by KaufmanC'^ ,,, 
aperture case are shown for comparison as solid symbols. The single 
aperture experimental results are observed to follow the same general trends 
as those predicted by the computer program solutions, but the experimental 
curves are shifted to lower perveance values. The magnitude of the pre- 
dicted minimum half angle agrees guite closely with that obtained experi- 
mentally. Unlike the theoretical solutions where this minimum occurs over 
a narrow perveance range, centered about a value of 1.5 x 10-s amp/volti/^ 

the Bxperimental results show a faiviu • • 

snow a fairly broad mimniuin, located around a mean 

perveance of 1.0 x 10 - amp/volt3/.. it is also noteworthy that theoretical 

solutions could be obtained at perveances up close to the Child's law limit 

maximum perveance that could be achieved experimentally by increasing 

thode em.ssion and hence thruster ion density was considerably less. 

figure 28 shows the effect of the variation in accelerator grid aperture 

diameter on the half angle vs. perveance plot for the single aperture case 

This figure shows that increasing the accelerator aperture diaumter to a 

value approaching the screen aperture diameter facilitates operation at 

higher perveance conditions. The one data point at a perveance beyond the 

theoretical Child's law limit of 3.03 x 10" amp/volt'/' is noteworthy 

This probably occurs because the plasma sheath position is „„t adeguately 

modelled by the parameter ^ u.;ed i„ tlie normalised perveance expre'.sion 
and defined in figure P3 
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Figure 29 shows tho effect of scroon grid thickness on the shapes of 
these curves is iniriiiiial for the single aperture case. The data presented 
in Figures 27 through 29 served the purposes of 1) verifying the adequacy 
of the computer ion optics solutions of reference [17] and 2) indicating 
what values of the geometrical variables would be of primary interest for 
the multi aperture experiments, 

Results obtained with mul tiaperture grid sets which more closely model 
actual thruster grid sets are presented as Figures 30 through 32 as a 
function of the same variables and parameters as those used for the single 
aperture grid data. Interactions between adjacent holes which occur when 
the multiaperture grid sets are used, appear to result in no pronounced 
departures from the general trends found with the single hole geometries. 
Here again a fairly broad perveance range at the minimum half angle was 
observed to be centered about a normalized perveance value of about 
1 X 10 amp/volt^/' . The half angle did not reach such high values at 
low perveances as it did in the single aperture cases. This effect may 
be due however to errors in i nterpretation of the ion beai 'files at 

low perveance where the ion beam profiles were generally difficult to 
analyze. 



Some general observations apparent from the multiple aperture work are: 



1. Decreasing the grid separation ratio results in a dramatic 
increase in the magnitude of the minimum half angle but does 
not change the shape of tlie curves si(|ni fi cantl y . 

2. Increasing the acc(? I orator t)nle (liametor at a constant ,•* , 

d^ 

I’c’sults in a signifiiant increase in the iiiaxiiiimn obtainable 
perveam.o iiut doc's not ctiaii'ie the minimum value' nl itu' hall atrili'. 
An ill I {’lerator aperture to 'a reen aperture (li.imeter i .iMo ol n,', 
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represents a limit below which the preveance capability of the 
grids is limited significantly. 

9 . d 

3. Reducing the screen thickness for a constant and -A , 

serves to broaden the minimum half angle region over a greater 
operating perveance range. An increase in the maximum obtainable 
perveance is also apparent, but a relatively small change is 
observed in the minimum half angles. 

t 

A point worth noting was that while a 30% reduction in from 0.185 

^s 

to 0.123 substantially broadened the minimum half angle range and increased 
the maximum obtainable perveance, only a slight improvement was noticed in 
both of these trends as the screen thickness ratio was reduced an additional 
50% to 0.062. This suggests reductions in screen thickness ratio below 
0.12 are not effective. The reason for this apparent lower limit is 
associated with the position and shape of the plasma sheath relative to 
the screen grid. For thick grids, ion recombination at the inner wall of 
the screen hole acts to reduce the ion density in this region, which dis- 
torts the sheath shape, reducing focusing. However it is possible that if 
the screen is too thin, the sheath can no longer attach itself to it and 
takes up a detached position slightly downstream of the grid. This slight 
reduction in the actual value of p,g, counteracts any focusing gains achieved 
by not having a distorted sheath edge. 

d. 

Plots of the beam diameter ratio ^ obtained with the single aperture 

^s 

configuration are presented in Figure 33. They exhibit correspondingly 
lower perveance values than those found in the computer solutions of 
Kaufman, although the general trends agree. There is a marked differ- 
ence, however, with the results for the multiple aperture grids shown in 
Figure 34. Unlike those of the single hole, the data spread for the three 
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FIGURE 34 
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different net-to-total accelerating voltage ratios used is minimal and 
the curves fall on essentially the same line. Also, the perveance values 
producing a given beam diameter are considerably lower for the multiaperture 
case. In both cases variation of the parameter ^ resulted in minimal 
scattering of the data. 

Conclusions 

f '1 

The computer solutions obtained by Kaufman for half angle variation 
agree qualitatively with experimental results. Minimum half angles pre- 
dicted by the computer solutions also agree with those observed experi- xi 

mentally although the perveances at which these minima occur do not agree. 

For all the grid geometries studied--both single and multiple hole— the 

minimum half angle occurs consistently around a normalized perveance of ] 

1 X 10"^ amp/volt3/2. The magnitude of the minimum half angle of an ion j 

i 

beam is reduced primarily as a result of increases in the grid separation j 

s 

i 

distance. The maximum obtainable preveance is increased primarily by i 

increases in the acceleration aperture diameter to values near that of the 

screen grid. Reduction in screen grid thickness tend to broaden the 

perveance range over which the beam half angle is a minimum down to the 

point where the screen thickness is about 0.12 times the screen aperture 

diameter. 

Charge exchange ion production within the ion beam makes absolute 

J 

ietermination of the beam half angle difficult. A better understanding of ; 

I 

the distribution of the charge exchange ion population needs to be ascertained i 

j 

from further experimental investigation in this area. I 


i 
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VISUAL OBSERVATIONS OF THE INTERIOR 
OF AN OPERATING HOLLOW CATHODE 

The hollow cathodes used in ion thruster applications consist of 
d tantalum tube capped by a thoriated tungsten plate perforated by a 
single small hole at the tube centerline. The mechanism by which this 
device supplies electrons to sustain an arc discharge when mercury is 
flowing through it is not well understood. Previous studies have shown 
however that a plasma does exist within the tube^^°^ and this suggests 
processes within the tube, which are not generally observable, play a 
significant roll in cathode operation. In order to facilitate observation 
of the interior of an operating cathode and hence observe phenomena that 
could lead to an understanding of cathode operations a 6.4 mm dia. hollow 

cathode was equipped with a quartz window and tested in a vacuum bell jar 
facility. 

Apparatus 

The apparatus used in this study and shown schematically in Figure 35 
consists of a 6.4 mm dia. hollow cathode supplied with an orifice plate 
having a 0.4 mm dia. orifice and fitted with a 1.6 mm thick by 1.27 cm 
diameter quartz window. Through this window the interior face of the 
orifice plate can be observed and photographed. The mercury flow rate 
through the cathode orifice was regulated by controlling the current to 
a standard SERT II vaporizer. The keeper was a wire loop positioned 1.5 inm 
from the cathode tip and the anode was a 3 cm diameter cylinder made from 
perforated stainless steel sheet metal. The cathode assembly was equipped 
with a tickler electrode to facilitate startup. The (athode ifcelf was 
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bare (no insert) and the interior surtaces were tree ot irregularities 
visible to the naked eye. The temperature at various locations on the 
interior face of the orifice plate was measured using an optical pyrometer. 

Emissive Material Deple tion 

Low work function oxides, generally in the form of Chemical R-500* 
are used inside a hollow cathode to facilitate startup and subsequent 
operation at acceptably low temperatures. There is evidence that this 
material may become depleted particularly if the cathode is operated at 
the high temperatures encountered at high arc current operating conditions. 
The phenomenological changes which accompany depletion are of interest and 
were observed in the following test. 

A new cathode with no insert was treated with a very small amount 
of R-500 (diluted tenfold before being added by a damp fibre dobber through 
the opening left when the window was removed). This cathode was then 
operated for about 135 hours at a high arc current (7A) and the interior 
orifice plate face was photographed and surveyed v/ith an optical pyrometer 
periodically to determine changes in its condition. Review of these data 
suggested the initial R-500 loading had been insufficient and had depleted 
after a few minutes of operation. A larger amount of R-500 was then added 
(undiluted from a damp fibre dobber). This addition caused a dramatic 
change in cathode operating parameters followed by a gradual return (over 
a 24 hour period) to those conditinns existing before the final R-500 
addition . 

The 'tfp(t of R-fiOO addition and depl('tion on ( athode operating 


* 
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paraiiiQters and cathode appearance; is shown in I inure' 36 and its accom- 
panyinc) captions. The captions show for example that tlie new cathode 
shown in the first column of photographs (probably without siqnificant 
R-500 but with thorium present in an undisturbed state) operates at an 
interior orifice plate temperature of 1370"C to 1200"C at an arc voltage 
of 17v and a keeper voltage of llv (the temperature range given is that 
observed from the edge of the orifice itself to the outer edge of the 
orifice plate). After the 135 hour operating period R-500 was added 
and the temperature dropped to 1190°C to lOOO^C, the arc and keeper 
voltages to llv and 4v respectively and the cathode appeared as shown 
in the second column. After another 24 hour operating period and R-500 
depletion the temperature rose to 1520°C to 1500°C, the arc and keeper 
voltages rose to 40v and 14v respectively and the cathode appeared as 
shown in the third column. These final conditions and the cathode 
appearance were similar to those observed just before the R-500 was 
added (after the 135 hour operating period). The photographs were ob- 
tained through a red #29 wratten filter and a 1" neutral density wratten 
filter using Ektachrome EX-135 film. Tite photographs can be interpreted 
if one considers the strong spectral littes of the elements present and 
the transmissivity of the filters at the wave lengths of these lines. 

The strong spectral lines listed in Table IV were determined by viewing 
the discharge through a hand spectroscope, and it was found that the 
depleted cathode showed only the merniry line, while the treated cathode 
shov/ed the additional line;. chatsK.fer i‘> Li c ol tlie Iwu'ium Unit i*. an active 
i n*)red i (; f 1 1 ot cliemical R-5('0. A review ol lahlc' IV 'aicpies Is thc' ted I iltet 
'diould paS) no i c|iii i i c .in t, I'.iclialion I rom tlie ineiiiiev. I he red pho i.oc|t dph'; 
should liowi'ver show Hie loi.iliun ol cs-s iled neiili.ils .md cit i nurse 
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thermal radiation from the hot cathode surface. The neutral density filter 
photographs should on tfie other liand stiow the location of the excited 
mercury neutrals whose r'adiation tends to dominate in the plasma. When 
some barium is present it can be seen on the cathode centerline (red 
filter photographs 1 and 2) but wlien it has been depleted the orifice 
appears black (red filter p!iotoc)raph 3). The neutral density photographs 
show the continuum r^adiation from tfie orifice plate becoming so intetise 
when the R-500 is deiileted (column ^) that it tends to overpower the 
mercury radiation and ’diow Uns)uuli as (n angc'. I he isrtteni ol tliis orange 
radic^tioti cati also bf* s('cn on (lie is ininal red 1 i 1 l('r photograpl)> i)ut it 
is not visible in ttn‘ reprodus 1 i ons includ(u;l lu'ce. Ibis patteni is 
a p[)ti tSMi 1 1 y related I e a va r i • t > i on 1 1 1 s in i a i ( ' 1 1 ■ s tor*’ n ^ I he in i i i < < ‘ p 1 a 1. 1 ' 
pf^odiicc'd dur i iri iii<uin \ i\< I m e . 

pho I ngi npii . irnr 0 iii'i slnii! .im- inc.rnl in lir' i-illiodi 
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region plasma when R-!iOO is present and this in turn would indicate barium 
migration takes place during cathode operation. An additional verification 
of barium migration was observed during the first few minutes of cathode 
operation after R-bOO addition. Initiallyj intensely radiating spots were 
observed through a red filter at locations where the R-500 was added. 

Some of these spots appear in the second red photograph of Figure 36 for 
example. After about 20 minutes of cathode operation these spots had 
disappeared and the orifice plate appeared to radiate quite uniformly 
through a red filter. 

Effect of Arc Current Variations 

At one point during the photographic study of an operating cathode the 
arc power supply current output was oscillating between 0.5A to 2A on a 
period of about a minute. Observation of the interior orifice plate sur- 
face through a red filter during this time showed that the entire orifice 
plate heated up when the arc current was high and then when the current 
suddenly dropped to one half ainpe'-e the radiation persisted for a period 
of about the thermal time constant before it decayed. The shape of the 
radiation patterns and the nature of their decay with time is suggested 
in the following sketch where early patterns (solid lines) and later 
patterns (dotted lines) are identified. This observation leads one to the 
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cont;lus1on that tho entiro inter ior oriliio plate partic:i[)ates in the 
eiiiisriion process ratlier than some small localized spot such as tlio orifice 
region. 

Cat hode Em ission Mec hanisms 

The Richardson-Dushman equation, which can be used to calculate 
thermionic electron emission from a surface, has generally yielded electron 
emission current densities based on measured external orifice plate 
temperatures that are too low to explain measured cathode emission currents. 
If one uses the internal orifice plate temperatures measured in this test 
series for an R-500 treated cathode (1100°C to 1200“C), a work function 
typical of a barium on tungsten surface (1.6 eV) and a Richardson equa- 
tion coefficient of 0.015 x 10'’ A/m-’"K’ the calculated emission current 
density is found to be 3.8 x 10“' to 11 x 10'' A/m^. Since the entire 
orifice plate was suggested in the preceding section as partici- 
pating in the emission process it is logical to use the interior surface 
area of the orifice plate as the emitting area and this leads to an emis- 
sion current of 1.2 to 3.4A. Since the cathode was operating at 7A this 
model again appears to be inadequate. A second experiment performed with 
an internal heated emitter provides however a clue to the source of addi- 
tional current. This test demonstrated that a heated emitter within a 
6.4 mm dia. cathode could effect an arc current to an external anode which 
was two to three times the emission cur'rent of the emitter when the emitter 
was held at cathode potential. Ttiis two to three fold amril i ficat ion ol ti)(' 
emission current sugcjests the calc ulatefl thf'nninni< emiv.ion <:urrc'nl would 
be amplified to an ar> (urreiil. level ol i’.-l (-o Id.fiA. ihi'. c a I ( ii I a1 ed 


current level does bracket the measured 7A current and this suofJf-*^>ts the 
emission mechanism in a hollow cathode involves thermionic emission coupled 
with a plasma amplification effect. As barium depletion occurs and the 
plate temperature rises to 1500°C to 16flO°C this model suggests the surface 
work function would have to increase to about 2 eV to effect the 7A arc 
current assuming the thoriated tungsten Richardson coefficient of 0.04 x 10^ 
A/m2(°K)2. This is reasonably close to the 2.7 eV work function of 
thoriated tungsten. 



CULU CATHOOII SIAKTIJP TFiSTS 


Cathodos whidi can bo started reliably from cold conditions should 
facilitate a reduction in the propellant los. during startup. A preliminary 
study of means of achieving mercury hollow cathode startup from tempera- 
tures below 100°C was therefore conducted. The methods used to approach 
this task are described and tests conducted with 6.35 and 3.2 mm cathodes 
are detailed in this section. Cathode starting performance using internal 
and external electrodes and flow variation was determined. Ranges of 
cathode starting parameters are discussed using these configurations. 

General App roaches to Cathode Starting 

There are several methods that are available for initiating cathode 
operation. Briefly, the basic methods are: 

1. Lower the work function: This approach is generally accomplished 

by materials selection. Low work function materials are provided on the 
surfaces of the cathode or impregnated in pores within the cathode com- 
ponents to insure the availability of this active material throughout 

the desired life of the cathode. 

2. Increase the voltage stress: There are several alternatives to 

this approach. Shaping and positioning of the electrodes can significantly 
increase the voltage stresses. Movable electrodes can also be used. 

Increasing the applied voltage to existing electrodes as well as intro- 
ducing high voltage pulses on auxiliary (or tickler) electrodes have 

been used. High freguency alternatin(> potentials can also be applied. 

3. Decrease the breakdown potential: !li()her cathode (low rates or 

mechanical ol.struction of t!i.- normal flow s.'rvi's to de< reaso the l>reak- 
down potential by incn-asing (he gas pn-ssuro in i he propellant flow field. 


P 
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4. Increase the temperature: For the purposes of these tests the 

cathode was to be cold, but it was assumed under this ground rule that 
small metal emitters that could be heated rapidly to emissive temperatures 
at low power levels would be acceptable. 

Apparatus and Procedure 

A 6,35 mm diameter hollow cathode having a 1 mm diameter orifice was 
modified to facilitate evaluation of several of the approaches described 
above. Figure 37 is a cutaway sketch of the first of the configurations 
tested. The emitter coil and spike axial anode were operated as a diode 
pair at various bias conditions and both electrodes were also operated and 
biased independently. Figure 38 is a cutaway of the second configuration 
tested. The axial, 0.25 mm diameter, 10 mm long tungsten emitter was 
operated using the co-axial cathode wall as the anode. As both of the 
figures suggest the cathodes were mounted on a fitting containing a quartz 
window to allow viewing the internal cathode surfaces and electrodes 
during testing. No insert was used but an R-500 coating was applied to 
the inside face of the orifice plate. An external keeper was used with 
both of these configurations as was a tickler electrode on which a 12 kv 
potential spike could be applied. These configurations allowed preliminary 
testing of a large number of startup concepts without altering the cathode 
emission surfaces as a result of exposure to air. 

The general test procedure required stabilization of the flow rate 
(vaporizer temperature) while keeping all cathode hardware above the con- 
densation temperature of mercury. This could be quite accurately accomp- 
lished by eliminating visible condensation from the internal cathode and 
window surfaces. The keeper potential was set at 40(1 volts and if an 
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internal emitter was used i I: was a.'t. ,ii. a powcir lovol tliat |)rovided one 
amp internal emission. c:aUiod(> Marl in. | irsl/. were then conducted, and 
a tentative best starting teclinif|ue was selected by making comparisons 
of the results. 

Final tests were then performed on the enclosed keeper 3.2 mm dia. 
cathode assembly shown in Figure 39. This configuration allowed testing of 
the tickler concept and an axially inserted internal electron emitter. 

Since the primary interest of the program was the 3.2 mm enclosed keeper 
cathode, tests conducted both the 6.35 mm and 3.2 mm dia. cathodes covered 
the 187 to 4 mA flow rate range considered typical of the 3.2 mm cathode. 

Results and Discussion 

The tests on the configuration .luivyri in figure 37 were brief and 
consistent. Biasing the internal giif.,- anode positive with respect to the 
cathode resulted in immediate vapnt izaiion of the anode. Using this anode 
as an internal tickler resulted in rai.id starts followed by the need for 
spike anode replacement. The emitter coil produced diode operation when 
heated and biased negative with respe. t to tite cathode wall and startups 
were achieved. Operation of tfie internal cmiitter was sporadic, however, 
with either very low cur’rent, fiigli vulta.).) (9UV) discliarcjes or heavy 
discharges at under 10 volts o cin i iiri. Ilu' (oil ,. of this emitter were 
very difficult to center witli i in liP .athnde wall and tiris is (,on- 

sidered to be the probable rau .r o: ini' imnii.i' leni nr.ults. 

The configuration sliown in I ignii' g.iv.' I lent starting .liaracter- 
i > t i (, ■> dut i n(| th(' initial I ms I .. ni i i > V/ . I ■ l. . | < mI is thr i n I erna 1 n » r.mgr- 
iiietit ior the primary serii". mI m i r' '• I n'l i n.n v i|•sl^ h.id stnpv/n Unit 


rel i<it)le '.tart'. . on M nsid i | v 
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emitter and the external tickler at cathode flow rates greater than 200 mA, 
and it was therefore decided that the most meaningful comparisons of these 
startup techniques would be obtained at flow rates below this value. 

The tests were begun by establishing the desired flow rate and then 
raising the cathode temperature above the propellant condensation limit. 

This limit was found to be at indicated temperatures ranging from 135 to 

«• i 

160°C. This band was inferred by observing droplet disappearance and 

vacuum system pressure surges which occurred when the vaporization temperature : 

was reached. Tests conducted with non-condensable propellents showed that 
room temperature ('v20°C) starts could be achieved routinely, but temperatures ; 

above the condensation limit were required before startup could be achieved | 

with mercury. i 

■I 

When the internal emitter was used for startup, emitter heater power I 

was switched on to the preset value to start the cathode. After the keeper | 

discharge lit (rarely longer than the time constant for emission from the ^ 

j 

internal cathode) the emitter heater power was shut off. At cathode flow 1 

rates less than 200 mA the cathode was self-extinguishing at keeper currents 
in the neighborhood of 0.2A without internal emission. After extinction, 

1 

startup was reinitiated immediately and these cycles were continued as the ’ 

tip temperature slowly rose from 150°C to 250°C. At 250”C starting cycles 
were terminated until the tip temperature had cooled to 150”C and the 
starting cycles were then repeated. These startup tests were conducted 

I 

for flow rates of I-/, 140, 95, 45, 30, 10 and 4 mA equivalent. At all of j 

I 

i 

the flow rates the cathode started every time the internal emitter was j 

turned Oi if the emission level was at least one ani[). ^ 

The effect of tip heater power’ was liriotly investigated at flow rates I 

below 10 mA. At ti|) temperature', above KIOO T opetsilion loulil be '.u .tained 
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occdsioiiJlIy lull. I.Im' .lis.lMrMo w.r, ikJ wil.lioul; iiiicrii.il oiniv.ion. 

Fi(jiiro 40A show*. I'lio v.it'i.iLinii nl kci'iipi' poicirl, i.i ! wiUi itil.cni.tl oini'.'. ion 
for several flow ral.(-. when I.Ik' | (mii i I: l.er was l.i.r.n'cl |() volts ne()a- 

tive. As the How rate was deerecn.e'l , the ititerriiil I'liiission required to 
maintain the keeper voltaqe at a(.i.e|it.at) |y low values increased. Each flow 
rate has a distinct miniiiiuiii keeper voltane (horizontal asymptote) at a 
given internal bias. Figure 40B shows the effect of internal emission on 
keeper voltage for 10 and 30 volts internal bias at the lowest flow rate 
tested (4 mA). The keeper voltage decrease with changing internal bias 
was greatest at the low flow rate and was iiardly noticeable at higher flow 
rates. In all cases, starting was critically dependent on sufficient 
internal emission. Although at iGav.t one aiii|) i'‘'ternal emission was re- 
quired to start the keeper disdiarge, operation of the catliode could be 
sustained at any flow rate tested at an internal emission of 150 inA. 

The capacity of the keeper power supply is, of course, critical to 
startup under many cathode conditions. The dOO volt statfing potential 
was low enough that cathode ignition would not oc(.ur reliably using the 
ticklei eleettode at flow ralnvi |iclni,v (dxjul. .10 mA. Ihi'i condition existed 
because there was not enoinih ener(|y ■;torc.| in the i.ircuit to sustain the 
discharge. Keeper potentials ol OOd to i:’ 0 (i volt, wc-re suKicient however 
to maintain operation ii.imi ,i Mdl.lcc sl,,nl al all ol the I low rates in- 
vest iga ted. 

fwo other i n trni's i. i ng clu'i i. . wino oP.mi veil during |||(i '.tarlitig tests. 
I) Oui ing hiijli tt.Miipeisi l urn oprr.il imi,, bar uim uppnngii ly iiiigiatr'd to (Im 

nun-opi.Mstl ini| inlnin,il nmiiir.|. liii v,m m Inl.ilrd I ibr n,.,\l nmillnr 

start alli'liipt v.'liidi Wn l nv n i l .P. i ■■ ni dn.i.ihi mu' Ini Wllh 1.1)11.11'/ 

dn. I !S|-. i I). I I'lM i . I nn .1 I .1 1 1 , : , . , . . , ■ 1 - I : 

' ’ * ' .id' ' . ' : * ! ' : I J I M ; I ‘ 1 I I t ( I I I ' / 
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A. AT lOV INTERNAL BIAS. 



B. AT 4 mA FLOW RATE 


EFFECTS OF INTERNAL EMISSION ON KEEPER VOLTAGE 


FIGURE AO 
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condensation limit, the intern.) I I'liii ti.ei’, due to its low mass and lessor 
thermal lao, collected iiiercioy. I his vois |K(stul.il.ed I roiii tlie next eiiiittei 
start attempt which would be an "instant on" start with a one decade 
pressure pulse of very short dmsil ion seen on tl)e vacuum instrumentation. 

A 3.2 mm diameter enclosed kec'per l atliode was also tested. The 
cathode was operated, as supplied, with no insert and tree of emissive 
material. At mercury flow rates of about 100 iiiA and keeper potentials 
of 500 volts, startup from cathode temperatures just above the condensation 
limit were demonstrated using a 0.25 mm diameter high voltage tickler. 

The tickler was positioned at the edge of the keeper orifice, 1 mm down- 
stream with the axis of the tickler wire parallel to the cathode axis. 

(See Figure 39). Attempts to achieve startup at lower flow rates or 
temperatures were generally not successful. Increasing the keeper potential 
above 500 volts resulted in breakdown thi'ougli the aluminum oxide keeper 
support tube. The 3.2 mm cathode su[)port elbow was modified to allow the 
introduction of an internal emitter is shown in Figure 39. The axial tube 
mounting method was required because of the small bore of this cathode. 
Because of the closely-spaced c|eoniet»y, the alumina emitter supports fused 
during each run. Operation could not be sustained for sufficient periods 
to allow data to be taken. Additional oi>eration with the 3.2 mm cathode and 
internal emitter was abandoned. Ilu' eng ineerinii difficulty associated with 
implementing the emitte)' in tliis .i/e ( liiiodi' w<is beyoixJ tlie scope of the 
present effort. 

Li fetjmc;_ 

A1 thouffli ('Xtendr'd I i 1 1 ' ti'.l . wi'M' nnl iim IimIimI in this ellnri, ftie 


cycling capabi 1 i t it' 
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to log the number of starts. After the preliminary tests were conducted 
and the procedures described previously were established, a new internal 
emitter was installed for the remainder of the tests. The emitter was 
used to make repeated cathode starts during the periods when the flow 
rates were being stabilized and measured. During the course of this pro- 
gram a single emitter accumulated 5100 starts. Failure of the emitter 
occurred during a Monday morning startup at a power level 1/4 to 1/3 of 
the normal operating power. None of the characteristics of the 0.25 mm 
diameter tungsten emitter had been observed to change prior to this failure. 
The primary significance of the cycle testing is that an unoptimized, 
first- try configuration achieved over 5000 cycles when operated in a vacuum 
system that was shut down every night over the ten day period of the tests. 

Conclusions 

Cathode starts can be achieved at temperature level: down to the 
condensation temperature of the mercury propellant prevailing at hollow 
cathode pressures (135 to 160°C). The most successful means of starting a 
6.35 mm diameter cathode was through the use of an internal 0.25 mm diameter, 
10 mm long tungsten emitter which was biased negative with respect to the 
cathode wall. Successful starting and operation was achieved at propellant 
flow rates down to 4 niA equivalent. 

An external high voltage (12 kV) tickler was the next best starting 
method for this cathode and the only method that worked successfully in the 
3.2 mm diameter cathode. Below a 30 mA flow rate limit, higher keeper 
potentials were required to sustain a tickler initiated start. 


ION THRUSTLK HOLLOW CAIIIODl Pl.ASMA CHARACTLRISTH;:. 

The plasma inside an operatiiuj 0 . 3 ? cm diaiiic'tnr hollow cathode has 
been probed by Fearn, Philip and In order to verily their 

results and also to determitie if the |)robe used in their study altered 
the characteristics of the plasma in their rather small cathode a similar 
test was conducted under this grant using a larger diameter cathode. A 
0.63 cm diameter cathode having a 0.38 mm diameter orifice was outfitted 
with a stationary Langmuir probe located on the cathode axis in the manner 
suggested by Figure 41. The cathode itself was supported by a structure 
provided by a quartz window in the manner suggested by Figure 35. The 
interior side of the cathode orifice plate was coated v/ith chemical R-500 
for the test, but no insert was installed within the cathode. 

The Langmuir probe was 0.13 mm in diameter and 0.46 cm long and was 
made of tungsten wire. It was biased using a battery circuit and the bias 
voltage and probe current (sensed across a lOOc; resistor) were fed to an 
X-Y recorder. Langmuir probe data were analyzed using least squares ex- 
ponential curve fits to the data in both the electron decelerating and 
accelerating voltage regions of the probe traces. A thin slieatii analysis 
was used in spite of the fact that the sheath thickness may be comparable 
to the probe diameter at some operating conditions. This assumption while 
it may introduce as much as two-1 old error in electron density should not 
affect electron temperature or plasma potential sicpii ficantly . 

Itesults 

As d.ita were (oll('cted and the tiiiK' t.tie proix' had been expo'. el to the 
plasma increased tlie pi-olu' ( harat t(>ri',ti' were ol)'.crve(i 1ri i hatiip' pro'.iim.ibl y 


TO BIAS AND SENSING 
CIRCUITRY 



0.63cm DIA. CATHODE 


CATHODE lANGMUIR PROBE 


FIGURE 41 
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because of probo contamination. Attompt'; to cloan tl)o probo by biasing 
it to either +45V or -4BV using a battory were not successful although 
the prooe was observed to be heated to incandescence by electron bombard- 
ment when it was biased +4BV. This contamination appeared to occur more 
capidly as arc current was increased and was presumably caused by con- 
stituents of the R-500 mix which was present. The contamination problems 
prevented exteried operation of the cathode and accurate determination of 
mercury flow rates but some data could be collected before contamination 
became seve ■ Figure 42 shows for example the variations in plasma 
properties with arc current determined from Langmuir probe traces before 
contamination became significant for an approximate mercury flow rate of 
100 mA. The magnitude of the electron densities and temperatures are 
observed to be in good agreement with those measured by Fearn et al.^^^^ 

The trends shown are considered qualitatively correct. 

Conclusions 

In view of the fact that essentially identical plasma properties have 
been measured in both 3.2 mm and 6.4 mm dia. cathodes it is concluded that 
the Debye snielding distance in an operating hollow cathode is sufficiently 
small to facilitate the collection of meaningful Langmuir probe data. Probe 
contamination due possibly to the constituents in the emissive iiii> present 
can introduce significant errors into the measurements. The plasma potential, 
electron temperature and electron density within a hollow cathode lie 
typically in the range of 4V, O.BeV and 10''Jc.nr' r('s[)ecti vely. 





HOLLOW (.AIIIODI Ml LIMD OKI I JU MIJDY 


It has been postu latdH i.lial I'UMhdn iMnissioi) cm (.ur*. a:, a rcy;ult 
of a surfacG GMiission iiUMHian i mm on tlif inl.oi ioi' ol an ion thruster 
hollow cathode operated on Muo'cury; tlniL L1 m‘ emission current density has 
a constant value in the ramje LU to ['>0 A/om- ; and that the emission occurs 
in the vicinity of tlie cathode orificeJ^"^^ As emissiofi current is in- 
creased then this mechanism simply requires an increase in the emitting 
surface area. Fearn, Philip and Pye^^^^-^ conducted an experiment utilizing 
a cathode orifice that was stepped to a larger diameter part way through 
the orifice plate. They found tlris sLo|.i introduced a corresponding transi- 
tion in the emission current/vol tage characteristic of the cathode. In 
order to verify this experiment and to observe visually the existence and 
motion of this emitting surface as arc current is increased a similar ex- 
periment was conducted. A cathode having tlie t'/O step orifice configuration 
shown in Figure 43 was constructed and install(^d on a fixture containing a 
quartz window in the manner suggestt'd in l ifuire 3'). Tlie cathode interior 
was not provided with a cathode insofM. but was coated with R500. 

Re sults 

Figure 44 preserits [)lot‘. cd ais < urnnit vs. arc voltage ol)tained in a 
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0.38 mm 
0.4 I mm 
0.53 mm 



STEPPED ORIFICE CATHODE 


FIGURE 43 
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ARC VOLTAGE ( V) 
A. HIGH FLOW RATE ( 250 mA ) 
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ARC VOLTAGE (V) 


B. INTERMEDIATE FLOW RATE (77 mA) 



ARC VOLTAGE (V) 

C. LOW FLOW RATE ( 23 mA ) 
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and current being determined by the current limiting characteristic of the 
arc power supply. At the intermediate flow (Figure 44B1 the arrow indi- 
cates the occurrence of an oscillation with a character again determined 
by the power supply. At the low flow rate (Figure 44C) the arc current 
continued to increase gradually with increases in arc voltage. The break 
in the curve of Figure 44B and the peak arc current condition of Figure 44C 
were characterized by a transition in the external discharge. At voltages 
above these values the plasma was observed to radiate throughout the bell 
jar whereas at lower voltages the plasma was observed only in the region 
between the anode and cathode. The data of Figure 44 show no evidence of 
any transition in the emission current voltage characteristics that might 
be associated with movement of an emitting surface to the geometrical steps 
in the orifice. 

Observation of the plasma within the hollow cathode through red and 
neutral density filters indicated the intensity of radiation from the 
plasma increased with arc current, but there was no evidence of an emitting 
surface that increased in area over the cathode orifice region as arc cur- 
rent was increased. 

Conclusio n 

Arc voltage/current characteristics and visual observations of the 
interior of an operating stepped orifice cathode do not support the theory 
of a constant current density emitting surface which exhibits an increase 
in area as arc current is increased. 
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APPENDIX A 


SINGLE AMD DOUBLE ION BEAM CURRENT DETERMINATION 
John R. Beattie 

Nomenclature 

Aq = grid viewing area (o =0), rn^ 

Ap = probe aperture area, m^ 

Ag = grid viewing area (o ^0), 

f(e) = correction factor 

i(r,e) = probe current, A 

I(r) = current crossing measurement plane, A 

j(r) = current density at grids, Am"^ 

J(R,o) = current density at measurement plane, Am"^ 

L = distance from grids to measurement plane, m 
r = radial coordinate measured from beam centerline, m 
r^ = grid radius, m 

R = radial coordinate measured from center of circular viewing 
area, m 

0 = probe pitch angle, 

= dispersion angle, ” 

“m=.v “ positive dispersion angle, ^ 

'min ” negative dispersion angle, ° 

!' = probe acceptance half-angle, ^ 
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I ntroductio n 

The determination of the single and double ion beam currents 
requires viewing a point on the grids from several different angles 
since the ions leaving the accelerator system follow trajectories which 
are not parallel. That is, the single and double ion dispersion profiles 
must be determined as a function of the thruster radius and this can be 
accomplished by the use of an articulating t x § momentum analyzer. 

This appendix presents the derivation of the equations which are used to 
transform the momentum analyzer output into integrated single and double 
ion beam current components. 

Beam Current Equations 

The t X momentum analyzer probe is used to measure the current 
leaving a small area of the grids as a function of the dispersion 
angle o. The probe is pointed at a fixed area of the grids at different 
angles and the current due to ions leaving this area at each of the angles 
is detected by the probe. When the probe is positioned at o = 0 the 
viewing area is circular as illustrated in Figure Al-B. With the probe 
at any other angle o the viewing area is elliptical as indicated in 
Figure Al-C. To account for the associated area change the current density 
at the analyzer plane is weighted by a factor f(o) which is equal to the 
ratio of the circular to elliptical areas. For simplicity in the deriva- 
tion which follows, the probe entrance aperture plane is assumed to 
remain a fixed distance L from ttie plane of the grids as illustrated in 
Figure Al-A. 

The ion (urrent density at the nijasureiiient [)lane •)(«,") is 
calculated I roiii l.lie pt ol>e currc'nt i(r, ) and fht! probe a|)erfure area A 
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by tho relationship 


j(R,,,) (Ai) 

P 

The current crossing the measurement plane as a result of ion injection 
from a circular region on the grids at a radius r is obtained by integrating 
this current density over the signal area at the measurement plane 


I(r) = 



where 


From Figure Al we see 



J(R.o) cos n dA 


dA - 2irRdR . 


R ^ L tan I) 


(A2) 


and 

dA = 2nL-' sec-’ o|tan oldo . (A3) 

Since the dispersion profiles are generally not symmetrical about o = 0, 
the integration is carried from *' = to o = and the result is 

divided by two to obtain the current I(r). 

The ion current density at the grids is obtained by dividing the 
current crossing the measuromctit plane by the area from which th(? cuMont 
left. That is, 

i(r) l(r)/A^^ (AT) 


ITie ion i un-f'iit tan lie oblaiiuMl by i n l.egra t, i n<| t tie . urrenl di'n'.ity ovoi 
Uie tp’itl a.rt'.i l>V ii'a' of fin' o /pn".'. i rm 
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Ig 2irrj{r)dr (A5) 

where is the grid radius. Implicit in Equation A5 is the assumption 
of azimuthal symmetry. Combining equations A1 , A2,. A3, A4, and A5 gives 
the integral equation for the beam current 


Ig = ri (r,e)f (o) |sin elsec^ e dedr 

— I 


(A6) 


P 0 


mm 


The function f(e) can be derived by use of Figures A1 and A2. 
The circular viewing area A^ is given by 


Aq = irL^- tan^ ijj . 


(A7) 


The elliptical area A. is 


irL^ tan \h [tan (e + lii) - tan (0 - ij 


2 cos j 

^sin“^ 

“ tan (e+ij;) + tan 

2 j 

} 


F.quation A8 can be simplified by use of the following approximate 


(A8) 


relation which introduces an error of less than 1% for dispersion angles 
less than about 25° 


1 irL^ tan i^[tan (o^lj;) - tan (o-iJj)! 
0 “ 2 cos 0 


(A9) 


The correction f(o) is obtained by combining Equations A7 and A9 to give 


f(n) = A„/A„ 

• 2 cos 0 tan t . (AlO) 

tan T'i+'iT tan (ii-'jr) 
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Data Acquisition a nd Red u ction T ec hni que 

Analyi'er probe current data i(r,o) were obtained in the following 
manner. The analyzer plate voltage was adjusted to allow detection of 
the desired ion current component, The probe pitch angle o was set to 
within 0.05° of the desired value with a stepper motor and the probe was 
swept vertically through the icn beam from r = -r^ to r = r^ by means of 
a motorized drive assembly. The output of a linear potentiometer was used 
to monitor the vertical position of the probe and this signal along with 
the analyzer output was recorded on an X - Y plotter. The probe pitch 
angle was then increased by a 5° increment and the procedure repeated. 
Since the beam is assumed axi symmetric, the probe angle can always be of 
the same sign and the positive or negative dispersion angles are de- 
termined by r > 0 or r < 0. This reduces the number of off-axis probe 
pitch angles by one half which reduces the data acquisition considerably. 

Data reduction consists of digitizing the plotted data for input to 
a digital computer program which performs the double integration of 
Equation A6 by numerical means first ^"o^ the doubly charged ion data and 
then for the singly charged ion data. The ratio of these two currents is 
then the double-to-single ion current ratio. Ion dispersion profiles are 
also determined using Equation A1 for the single and double ion data. 


Af’PENDIX B 


LANGMUIR PROBE ANAEYSIS 
John R. Beattie 

The Langmuir probe data presented in this report were analyzed by 
the numerical procedure of Reference [10]. These results generally 
indicate higher Maxwellian electron temperatures and lower primary to 
Maxwellian electron density ratios than earlier results which were ob- 
tained by graphical methods. The purpose of this appendix is to 
explain these discrepancies and to point out the large errors the 
graphical procedure can introduce when the electron temperature is 
greater than about 4eV. 

In the retarding field portion of the Langmuir probe trace the 
two-group plasma assumption predicts the following current-voltage 
relationship^^*^^ 


I -- Bi + BvV + B;j exp (B 4 V) (Bl) 

where I is the probe electron current, V the probe voltage, and the 
constants B^. are determined by the plasma properties. Introducing 
these properties, the expression above can be rewritten in the following 
nondimensional form 


I 


r-=^ 

sat ^m 


'4T 

m p 


(s 


+ V) + exp 


V - •!■ 




(BZ) 


where n^^ and n^^^ are the primary and Maxwellian electron densities, 
t) ^in primary electron energy and Maxwellian electron tempo 


:ira- 


ture, ! the plasma potential, and the Maxwellian electr 


on 


Stiturdtion current. l^r}()rn'.fMit dt i vn viluf“. ni tlir propnr't i'"' ",/ti , 


... , and <|.p are presented in Table B1 for a mercury electron-bombardment 
thruster. Usinq these values. Equation B2 was used to calculate 

TABLE Bl. ASSUMED PLASMA PROPERTIES 


Property 

Assumed Value 


5% 


30eV 

- 

40V 

.L— 


dimensionless current as a function of probe voltage for various electron 
temperatures. The results are presented in Figure Bl which indicates the 
variation of total and primary electron current with probe voltage for 
three values of the Maxwellian electron temperature. Figure Bl-A corres- 
ponds to a fairly low electron temperature of 2eV and we see the primary 
and total current lines are tangent (i.e., the Maxwellian current is very 
small) over a large voltage range. Graphical analysis of these data 
should therefore yield good results. Figure Bl-B is for an electron 
temperature of 6eV and we see at this higher temperature the two curves 
are never tangent. Analysis of these data by the graphical technique would 
result in erroneous plasma properties since the method requires the primary 
current line to be tangent to the total current line. To estimate the 
amount of error one might expect, the dashed curve labeled A-A in Figure 
Bl-B was drawn tangent to the total current line in such a manner as to 
(jive the correct primary electron energy. Analysis of these data using 
the line A-A as the |>riiiiary current re-.ulted in errors in the primary 
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electron number density, Maxwellian electron temperature and density of 
160';.'., 20'Z, and 1',:'., respectively. Figure Bl-C indicates how the problem 
becomes even more severe when the electron temperature is increased to 
10 eV. The line labeled A-A in this figure is drawn tangent to the total 
current curve and gives the correct primary electron energy. However, in 
the accelerating field region of the trace the slope of line A-A is greater 
than the slope of the total current curve — a condition which is not reason- 
able physically. In order to obtain a realistic solution in the accelerating 
field region and still satisfy the tangency condition, a primary current 
line such a B-B might be used but this also introduces an error in the 
primary electron energy. Graphical analysis of Figure Bl-C using the 
line B-B as the primary electron current resulted in errors in the primary 
electron energy and density, and Maxwellian electron temperature and density 
of 13S, 307%, 35%, and 9%, respectively. 

The numerical procedure for analyzing Langmuir probe data does 
not impose the tangency condition of the graphical method and, as a result, 
does not introduce the large errors described above. For comparison pur- 
poses the total current curves of Figure B1 were digitized and analyzed by 
the numerical procedure. The plasma properties determined in this manner 
were found to be in error by less than 15' in the worst case. 

The curves of Figure B1 indicate the graphical procedure will intro- 
duce significant errors in the plasma properties when the Maxwellian 
electron temperature is greater than about 4eV. Ihe error introduced will 
in general result in iiigiier primary electron densities ami lower Maxv/el 1 i an 
elec.tron teiii|M'ra lures Uian l.lie true valuf* .. 




.4 



REFERENCES 


1. Kerslake, W. R., et a1 , "SERT II: Mission, Thruster Performance, 

and In-Flight Thrust Measurements," Journal of Spacecraft and Roc kets, 
Vol . 8, No. 3, March 1971, pp 213-224. 

2. Vahrenkamp, R. P., "Measurement of Double Charged Ions in the Beam 
of a 30-cm Mercury Bombardment Thruster," AIAA Paper No. 73-1057, 
October 31 - November 2, 1973. 

3. Peters, R. R. and P. J. Wilbur, "Double Ion Production in Mercury 
Thrusters," AIAA Paper 75-398, March 19-21, 1975. 

4. Kieffer, L. J., "Electron Impact Ionization Cross Section Data for 
Atoms, Atomic Ions, and Diatomic Molecules," Rev. Mod. Physics , 

Vol. 38, No. 1, pp. 15 -23, 1966. 

5. Shpenik, 0. B. and Azpesochnyi , I. P., "Excitation Cross Sections 
near the Threshold for Electron Atom Collisions," Optics and Spect ., 
Vol. 23, pp. 7-10, 1967. 

6. Kupriyanov, S. E. and Z. Z. Latypov, "Ionization of Positive Ions by 
Electrons," Soviet Physics J ETP , Vol. 19, No. 3, pp. 558-559, Sept. 
1964. 

7. McConnel, J. C. and B. L. Moisewitsch, "Excitation of Mercury by 
Electrons," J . Phys . B . , Vol. 1, No. 3, pp. 409-412, 1968. 

8. Gryzinski, Michal, "Classical Theory of Atomic Collisions. I. Theory 
of Inelastic Collision," Ph ys. Rev. , Vol. 138, No. 2A, p. A341 , 

April 19, 1965. 

9. Wilbur, P. J., "An Experimental Investigation of a Hollow Cathode 
Discharge," NASA CR-120847, December, 1971. 

10. Beattie, J. R., "Numerical Procedure for Analyzing Langmuir Probe 
Data," A IAA Journal , Vol. 13, No. 7, July 1975, pp. 950-952. 

11. Wilbur, P. J., "15 cm Diameter Ion Thruster Research," NASA CR-134755, 
December, 1974. 

12. Beattie, J. R., and P. J. Wilbur, "15 cm Cusped Magnetic Field Mercury 
Ion Thruster Research," AIAA Paper No. 75-429, March 19-21, 1975. 

13. Knauer, W., R. L. Poeschel , and J. W. Ward, "Radial Field Kauliiian 
Thruster," ^urnal o/ Spacec 2 'aJ\ _and__Roci;_ets , Vol. 7, No. 3, 

March 1970, pp.' 2r4T(-'25l . 

14. Kaufman, H. H., "1 xi'cr imental Invostinal ions of Arfioti and Xenon Ion 
Sources," NASA CR- 1 34;’,4S , June, 1 )/S. 




W -W- I ' n » m7 i ^ ■ f'T w iTy? 1 gtf • 






1^-, Wilbur P j., "txperimental Investiqation of a Throttlable lb an 
^ ' Hollow’cathod; Ion Thruster," NASA CR-1P1038, December, 1972. 

16 Kaufman, Harold R., "Ion-Thruster Propellant Utilization," 

NASA TN D-6591, December, 1971. 

17 Kaufman, Harold R., ■'Accelerator-System Solutions for Electron-Bombard- 
ment Ion Sources,” AIAA Paper 75-430, March 19-21, 1975. 

18. Masek, T. 0., "Plasma 

Thrusters," AIAA Paper No. 69-256, March J-b, lyoy. 

19. Kaufman, Harold, "Charge-Exchange Plasma Generated by an Ion Thruster, 
NASA CR-134844, June 1975. 

20. Fearn, D. G., C, «. PM'i^o^ate^s'lo; S^nTer°cS^lon^^ril^ers." 

22/23 Juni 1971, Braunschweig, West Germany, pp. 61-83. 

21. Anon., "Handbook of Chemistry and Physics," 36th Edition, pp. 2654-2658, 
Chemical Rubber Publishing Co., 1954-1955. 

22. Anon., "Handbook of Chemistry and Physics " 36th Edition, pp. 2722-2733, 
Chemical Rubber Publishing Co., 1954-1955. 

. . ,, n tinri 1^ I Geiler. "Probe Measurements of Discharge 

?n"iiS9 n;ct?onU.ba?Lnt’Enpine," Vol . 1 , 

No. 8, August, 1963, pp. 1815-1823. 


